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Abstract 

There is evidence for temperature fluctuations in Planetary Nebulae and in Galactic 
Hii regions. If such fluctuations occur in the low-metallicity, extragalactic Hii regions 
used to probe the primordial helium abundance, the derived ^He mass fraction, Yp, 
could be systematically different from the true primordial value. For cooler, mainly high- 
metallicity Hii regions the derived helium abundance may be nearly unchanged but the 
oxygen abundance could have been seriously underestimated. For hotter, mainly low- 
metallicity Hii regions the oxygen abundance is likely accurate but the helium abundance 
could be underestimated. The net effect is to tilt the Y vs. Z relation, making it flatter 
and resulting in a higher inferred Yp. Although this effect could be large, there are 
no data which allow us to estimate the size of the temperature fluctuations for the 
extragalactic Hii regions. Therefore, we have explored this effect via Monte Carlos in 
which the abundances derived from a fiducial data set are modified by AT chosen from 
a distribution with < AT < AT^ax where ATmax is varied from 500K to 4000K. It is 
interesting that although this effect shifts the locations of the H ii regions in Y vs. O /H 
plane, it does not introduce any significant additional dispersion. 
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1 Introduction 



The primordial abundance of '^He is key to testing the consistency of the standard hot 
big bang model of cosmology and to using primordial nucleosynthesis as a probe of par- 
ticle physics (Steigman, Schramm & Gunn 1977). The availability of large numbers of 
carefully observed, low-metallicity H ii regions has permitted estimates of the primordial 
helium mass fraction, Ip, whose statistical uncertainties are very small, ~ 1% (see, Olive 
& Steigman 1995 (OS), Olive, Skillman & Steigman 1997 (OSS) and references therein). 
However, there remains the possibility that in the process of using the observational data 
to derive the abundances, contamination by unacknowledged systematic uncertainties 
has biased the inferred value of Ip. Observers have identified many potential sources of 
systematic uncertainties (Davidson & Kinman 1985; Pagel et al. 1992 (PSTE); Skillman 
et al. 1994; Izotov, Thuan k Lipovetsky 1994, 1997 (ITL); Peimbert 1996) and, where 
possible, have designed their observing programs to minimize such uncertainties and/or 
to account for them. It is expected that the contributions from many of the poten- 
tial sources of systematic uncertainty (if present at all) would vary from Hii region to 
Hii region and from observer (telescope/detector combination) to observer, introducing 
along with a systematic offset in the derived value of Ip, an accompanying dispersion in 
the helium data. More insidious would be systematic errors in, for example, the atomic 
physics used to convert the observed equivalent widths to abundances, since an offset 
from the "true" value of Ip would not be accompanied by an enhanced dispersion in the 
fits to the data (e.g., Y vs. 0/H or the weighted mean of Y for the lowest metallicity 
Hii regions; see OSS). In this paper we explore one potential source of systematic un- 
certainty in deriving abundances from emission-line data for extragalactic, Hii regions: 
temperature fiuctuations. 

The empirical method to derive chemical abundances in an emitting gas has been 
used widely since it was first proposed by Peimbert and Costero (1969). The electron 
density and the gas temperature T are obtained from emission-line intensity ratios and 
are used to calculate the line-emissivity which, along with the observed line-intensities, 
provide the fractional abundance of the ion and, subsequently, the empirical abundance 
of a given element (see, for example, Osterbrock 1989). In low density regions, such 
as those in Hii galaxies, the permitted-line emissivities (like those of H and He) de- 
crease (slowly) with T and are usually independent of Ug, unless coUisional effects are 
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important. On the other hand, the forbidden-hne emissivities (hke those of O, N, S) are 
strongly dependent on T, and the dependence on the density may also be important. 
In Hii galaxies, the [S II] hne ratio indicates low electron densities (ne < 500 cm~^), 
and except for those with high temperature, the coUisional effects are negligible. Thus, 
the derived heavy element abundances depend mainly on a good determination of the 
temperature. For these regions, the temperature used for the high-ionization lines is 
that obtained from the [O III] hne ratio, Tqiii, while for the low-ionization hues the 
temperature is derived from Tqiii using results from photoionization models (see for 
example, PSTE). 

A similar method is usually applied to planetary nebulae. In contrast to H ii regions, 
for several PNe the Balmer temperature (obtained from the observed Balmer discon- 
tinuity) is also determined. In many cases is found to be lower than Tqiii (Peimbert 
1971; Liu & Danziger 1993). As pointed out by Peimbert (1971), this discrepancy could 
be due to temperature fluctuations, which, however, are not reproduced by the stan- 
dard photoionization models for these nebulae (Liu & Danziger 1993). If, indeed, the 
gas temperature is overestimated by Tqiii, with the true temperature given by T^, the 
heavy element abundances derived from forbidden lines using Tom are underestimated 
(Viegas & Clegg 1994). Indeed, the abnormal chemical abundances inferred for some 
Hii galaxies with WR features may result from just such an overestimate of the gas 
temperature (Esteban & Peimbert 1995). A well-observed giant Hii region, NGC 2363, 
shows just such a discrepancy between the Paschen temperature and Tqih in two knots 
( Gonzalez- Delgado et al. 1994). For the Hii galaxies used to derive the primordial 
helium abundance, such an underestimate of the true oxygen abundance may introduce 
a systematic bias into the inferred value of Yp. Furthermore, temperature fluctuations 
also will have a direct effect on the helium abundance determined from the recombina- 
tion hues (Esteban & Peimbert 1995). Although the ratio of He and H emissivities is 
very nearly independent of temperature, the helium-line intensities are modified from 
the pure recombination case by coUisional excitation (Cox &: Daltabuit 1971; Ferland 
1986) which is temperature-sensitive. If the "true" Hii region temperature has been 
overestimated, the coUisional correction has been overestimated and the "true" heUum 
abundance underestimated. This effect will be largest for the hottest, metal-poor Hii re- 
gions. Thus the combined effects, which tend to increase the derived oxygen abundance 
in the metal-rich H ii regions and the helium abundance in the metal-poor regions (see 
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Figure 1), will "tilt" the inferred Y versus 0/H relation, flattening the slope and increas- 
ing the intercept, Yp. In this paper we analyse these effects in an attempt to quantify 
the corresponding systematic uncertainty in Yp. 

2 The Method 

If the Paschen or Balmer temperatures were available, the corrections described above 
would be straightforward to implement. Unfortunately, for the low-met allicity, extra- 
galactic H II regions used to derive the primordial abundance of helium, there appears 
to be no data on the Balmer or Paschen temperatures. To estimate the magnitude of 
the possible corrections we have, therefore, adopted a Monte Carlo approach to quanti- 
fying the potential systematic uncertainty introduced into the determination of Yp by 
uncertainty in the temperature of the Hii region gas. 

We begin by adopting a "fiducial" data set using 43 Hii galaxies from PSTE, ITL, 
Skillman & Kennicutt (1993) and Skillman et al. (1994). The role of this fiducial set is 
merely to provide a comparison, in order to quantify the changes induced by temperature 
differences. For this reason we have not employed the unpublished data of Skillman et 
al. (in preparation, see OS and OSS) nor those Hii regions listed in PSTE which they, 
themselves, did not observe. We follow ITL and eliminate those regions for which their 
data is noisy or otherwise uncertain. For each Hii region in our fiducial data set we 
adopt the same algorithms used by these authors to derive the temperature and the 
fractional abundances in the low-ionization zones. The electron density is obtained from 
the [S II] line ratio and we use the Brocklehurst (1972) recombination coefficients for the 
He and H lines along with the collisional correction for the He lines from Clegg (1987). 
The fractional abundance of He+ is the (unweighted) average of the abundances obtained 
from Hel A4471, Hcl A5876 and Hel A6678. Ignoring any statistical uncertainties, we use 
these fiducial abundances to find a "standard" linear fit to Y versus 0/H (unweighted) 
in the absence of temperature fluctuations. We emphasize that we are not so much 
interested in the "best fit" Y versus 0/H relation for our fiducial data set but, rather, 
the differences in the Y versus 0/H relations between AT zero and nonzero (AT = 
To///— T). In Figure 1 each point in our data set is shown in the Y versus 0/H plane 
in the absence of temperature fiuctuations (AT = 0) and for AT = 2000K, joined by a 
solid line. As anticipated, the higher metallicity regions tend to be cooler (see Figure 
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2). This has the effect that the derived oxygen abundance in high-metalhcity regions is 
very sensitive to temperature ffuctuations as is seen clearly in Figure 1. 

For planetary nebulae the observed temperature difference AT = To///— T^ is typ- 
ically 2000K (Liu & Danziger 1993) with one PN having AT = 6000K. For the giant 
Hii region NGC 2363, AT > 3000K (Gonzalez-Delgado et al. 1994). Thus, we have 
recalculated the helium and oxygen abundances for each Hii region with AT chosen 
from a distribution which ranges from zero to AT^a^;, according to a probability P(AT) 
to be described below; AT^^ax is varied in the range 500 < ATj^ax < 4000K. We repeat 
this procedure 10,000 times for each choice of ATmax- For most Hii regions Y does not 
decrease with decreasing temperature as might have been expected from the effective 
recombination coefficients alone (Esteban & Peimbert 1995) provided we account for the 
coUisional effect on the He hues; for Hii galaxies with high gas temperature, the coUi- 
sional effect dominates, leading to an increase of Y. As seen in Figure 1, the decrease in 
T leads to an increase in the oxygen abundance. For each of the 10,000 reahzations (for 
each choice of ATj^ax) of our H ii region data set we fit a hnear Y versus 0/H relation to 
derive Yp and the slope AY/AZ, where the heavy element mass fraction Z 20(O/H) 
(see PSTE). 

3 Results: Y Versus Z 
3.1 Full Data Set 

The results for our full data set of 43 Hii regions, selected to have 0/H < 1.5x 10~^, are 
shown in Figures 3-5. As anticipated, the intercept, Yp, is systematically increased by 
an amount which depends on AT^a^; (see Figures 3a & 5). This increase in the intercept 
of the Y versus Z relation is accompanied by a flattening of the slope (see Figures 3b & 
5). To probe the sensitivity of our results to the adopted probabihty distribution of AT 
values, P(AT), we have considered three, simple distributions: flat, linearly increasing, 
linearly decreasing. The results for these three, different distributions are shown for 
AT^aa; — 4000K in Figure 4. Clearly, the effects are closely similar for all choices; in our 
subsequent discussion (and, in Figures 3 & 5) we present results for the flat distribution 
(P(AT) independent of AT). 

In Figure 5 is shown how the magnitude of the systematic offsets in slope and in- 
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tercept scale with AT^^aa;. For the not unreasonable choice of AT^a^; = 2000K, the 
increase in the inferred value of Yp is significant, comparable to some estimates of the 
upper bound to the systematic uncertainty in Y (OS; OSS). For larger values of AT^a^;, 
the systematic offset will be even larger. It might have been anticipated (see OSS) that 
such large systematic offsets would be accompanied by an increase in the dispersion 
of the abundances around the best fit Y versus Z relation. To test for this, for each 
realization we have computed cr, the variance of the residuals between the data (Y) and 
those values predicted by the corresponding fit for that realization {Yfn). However, as 
seen in Figure 5, for temperature fluctuations this effect of added dispersion, if present 
at all, is very small compared, for example, to the typical uncertainty in the individual 
Hii region Y determinations which are of order 0.010 (OSS). Thus, until there are obser- 
vations of H II region temperatures determined from hydrogen lines, it is difficult to set 
a firm upper bound to the magnitude of the systematic offset in Yp due to the possible 
presence of temperature fluctuations. It is, therefore, important to consider how best 
to analyze current data so as to minimize the potential importance of such temperature 
fluctuations. 

3.2 Truncated Data Set 

A signiflcant contribution to the systematic offset in Yp in the presence of temperature 
fluctuations comes from the flattening of the slope of the Y versus Z relation due to the 
large increase in oxygen abundance for the cooler, higher-metallicity Hii regions (see 
Figures 1 & 2). Therefore, the uncertainty in Yp derived from a flt to the trend of 
Y with 0/H might be minimized if the data set is restricted to the very lowest metal 
abundance Hii regions which are hotter. To explore this, for each realization of our 
Monte Carlos we identify the subset of Hii regions with low-metaUicity: 0/H < 0.9 x 
10~^. For these subsets we flt linear Y versus Z relations and compare the slopes and 
intercepts (as well as a) to those found for the corresponding low-metallicity subset 
in the absence of temperature fluctuations. The results are shown in Figures 6 & 7. 
As expected, the changes in slope and intercept for this low-metallicity set are much 
smaller. Indeed, the trend is even opposite that for our full data set: lower intercept, 
higher slope (although the effect is so small as to be only marginally signiflcant). The 
lesson is clear. If we wish to minimize the uncertainty due to temperature fluctuations 
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in Yp derived from a fit of Y versus Z, we should restrict our attention to the most 
metal-poor H ii regions. 

3.3 Lowest-Y Hii Regions 

Since it is generally accepted that the helium abundance has only increased since the 
big bang, an alternate approach to using Hii regions observations to infer the primor- 
dial abundance is to compute the mean of Y for those regions with the lowest helium 
abundances; then, Yp < <Y> (OS; OSS). This estimator is hkely to be robust against 
the systematic shifts due to temperature fiuctuations. To test this we have first chosen 
the ten lowest-Y H ii regions in our fiducial data set (AT = 0) and found the unweighted 
average of Y, <Y>io(0). Then, from our Monte Carlos, for various choices of ^Tmax we 
determine (for each realization) the ten lowest-Y H ii regions (often, but not always, the 
same as in the fiducial set) and compute A<Y>io = <Y>io(AT) — <Y>io(0). The 
distributions of A<Y>io for lS.Tmax = 2000K and 4000K are shown in Figure 8 and 
the trend of A<Y>io with AT^aa: is shown in Figure 9 where the corresponding change 
in the dispersion around the mean {a) is also shown. As expected, although temper- 
ature fiuctuations tend to increase <Y>io systematically (mainly due to the reduced 
correction for coUisional excitation), the effect is quite small. 

4 Discussion 

Temperature fiuctuations in the low-metallicity, extragalactic Hii regions used to in- 
fer the primordial helium abundance will lead to systematic offsets in the helium and 
oxygen abundances derived for those regions. In the absence of direct observations of 
the hydrogen temperatures for these regions we have attempted to quantify the effect 
of such fiuctuations on the derived value of Yp by performing a series of Monte Carlo 
calculations where the temperature fiuctuation for each H ii region is chosen from a dis- 
tribution whose maximum value is varied. Yp determined from a linear fit to Y versus 
0/H is especially susceptible to this source of error which fiattens the Y versus Z relation 
resulting in a higher intercept. This potential systematic error is insidious in the sense 
that the offset in Yp may be significant without a noticeable increase in the dispersion of 
the individual H ii region Y- values around the best fit Y versus Z relation. To illustrate 
the possible offsets which may be consistent with extant data, consider the results of 
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OSS for their full data set of 62 Hii regions with 0/H < 1.5 x 10"^: Yp = 0.234±0.005 
(0.005 is a 2-(T estimate of the statistical uncertainty). For AT^aa; = 4000K we should 
increase this value by 0.008±0.004 (see Figure 5). Adding the (uncorrelated) systematic 
and statistical uncertainties in quadrature will lead us to Yp = 0.242±0.006 (95%CL). 
Such a large correction to Yp would amehorate the "tension" in BBN (Hata et al. 1995) 
between primordial helium and the low abundance of deuterium suggested by Galactic 
observations (see, e.g., Dearborn, Steigman & Tosi 1995 and references therein) and 
by some direct detections of deuterium in low-metallicity, high-redshift QSO absorbers 
(lytler. Fan & Buries 1996; Buries & Tytler 1996). 

Until data on temperature fluctuations in these key extragalactic H ii regions become 
available, the best strategy is to analyze the current data in a manner which minimizes 
this potential systematic error. One possibility is to restrict attention to the lowest- 
metaUicity regions available. OS and OSS have done this for the subset with 0/H < 
0.9 xlO-^ for which OSS derive from a hnear fit to Y versus 0/H: Yp = 0.230±0.007 
(95%CL). From our Monte Carlos we find for this subset that the systematic correc- 
tion is hkely very small (indeed, it may even be negative!); for AT^a^; = 4000K, AYp 
= — 0.004±0.005 (see Figure 7). If, perhaps naively, we apply this correction to the 
OSS result, we infer: Yp = 0.226±0.009 (95%CL). In this case the "tension" between 
primordial helium and low primordial deuterium is exacerbated. 

Another "safe" approach to using the existing data to derive an estimate of Yp is 
to ignore the metallicity information and simply take a mean (<Y>) of the helium 
abundance for the lowest abundance Hii regions. Since helium is only expected to 
increase after BBN, this provides an upper bound to Yp. In general, the Hii regions 
with the lowest values of Y tend to be the lowest-metaUicity regions which are also the 
hottest (see Figure 2) . For such regions there has often been a non-neghgible correction 
for collisions in deriving Y from the helium line intensities. If the gas is, in fact, cooler, 
this correction has been overestimated and the "true" value of Y should be larger. Thus, 
temperature fluctuations will increase <Y>. From our Monte Carlos we have selected, 
for each realization, the ten H ii regions with the lowest Y- values and we have found, for 
AT^aa; = 4000K, a small systematic increase: A<Y>io = 0.002±0.001 (see Figures 8 
& 9). For their ten lowest-Y regions, OSS find <Y>io = 0.230±0.006 (95%CL), so that 
even with our largest correction we infer a revised mean of 0.232±0.006 suggesting that 
Yp < 0.238 (95%CL). Here, too, the "tension" between D and ^He fails to be relieved. 
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5 Summary 



Temperature fluctuations in low-metallicity, extragalactic H ii regions may have a signif- 
icant effect on the determination of the primordial abundance of heUum. If present, they 
may increase the metalUcity of the higher-metalhcity, relatively cooler regions and in- 
crease the helium abundance of the more metal-poor, hotter regions, tilting the inferred 
Y versus Z relation and leading to a higher, zero-metallicity intercept (Yp). Such a sys- 
tematic offset is not necessarily accompanied by a signiflcant increase in the dispersion 
of the data around the best flt linear Y versus Z relation and, therefore, may remain 
invisible in the absence of direct data on temperature differences in these regions. It 
is clear that such data is crucial to constraining the uncertainty in Yp. In the absence 
of such data, we have noted that restricting attention to the lowest-metallicity regions 
(Z < Zq/10) will tend to minimize this systematic offset. Alternatively, the mean of 
the helium abundances of the lowest- Y Hii regions is also robust against the effect of 
temperature fluctuations. 

Acknowledgments 

S.M.V. and R.G. owe a debt of gratitude to P. Benevides-Soares for valuable discus- 
sions. The work of S.M.V. and R.G. is partially supported by FAPESP and CNPq in 
Brasil; that of G.S. is supported by the DOE. This work was initiated when S.M.V. and 
G.S. were visiting the INT at the University of Washington and they are grateful for the 
hospitality and stimulating environment provided there. 

References 

Brocklehurst, M. 1972, MNRAS, 157, 211 
Buries, S. & Tytler, D. 1996, ApJ, 460, 584 
Clegg, R.E.S. 1987, MNRAS, 229, 31P 
Cox, D.P. & Daltabuit, E. 1971, ApJ, 167, 257 
Davidson, K. & Kinman, T.D. 1985, ApJS, 58, 321 
Esteban, C. & Peimbert, M. 1995, A&A, 300, 78 
Ferland, G.J. 1986, ApJ, 310, L67 



9 



Gonzalez-Delgado, R.M., Perez, E., Tenorio-Tagle, G., Vilchez, J. M., Terlevich, E., 
Terlevich, R., Telles, E., Rodriguez-Espinosa, J. M., Mas-Hesse, M., Garcia- Vargas, 
M. L., Diaz, A. I., Cepa, J. & Castaneda , J. 1994, ApJ 437, 239 

Hata, N., Scherrer, R.J., Steigman, G., Thomas, D., Walker, T.R, Bludman, S. & 
Langacker, P. 1995, Phys. Rev. Lett., 75, 3977 

Izotov, Y.L, Thuan, T.X., & Lipovetsky, V.A. 1994 ApJ 435, 647. (ITL) 

Izotov, Y.L, Thuan, T.X., & Lipovetsky, V.A. 1997, ApJS, 108, 1 (ITL) 

Liu, X. & Danziger, J. 1993, MNRAS, 263, 256 

Ohve, K.A., Skillman, E. & Steigman, G. 1997, ApJ, In Press (OSS) 

Olive, K.A., & Steigman, G. 1995, ApJS, 97, 49 (OS) 

Osterbrock, D.E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei, 

University Science Books, Mill Valley, Cahfornia 
Pagel, B E.J., Simonson, E.A., Terlevich, R.J. & Edmunds, M. 1992, MNRAS, 255, 325 

(PSTE) 

Peimbert, M. 1971, Bol. Obs. Tonantzintla y Tacubaya, 6, 29 
Peimbert, M. 1996, Rev. Mex. Astr. Astrofis., Serie de Conferencias, 4, 55 
Peimbert, M. & Costcro, R. 1969, Bol. Obs. Tonantzintla y Tacubaya, 5, 3 
Skillman, E., & Kennicutt, R.C. 1993, ApJ, 411, 655 

Skillman, E., Terlevich, R.J., Kennicutt, R.C., Garnett, D.R., & Terlevich, E. 1994, 
ApJ, 431, 172 

Steigman, G., Schramm, D.N. & Gunn, J. 1977, Phys. Lett, B66, 202 
lytler, D., Fan, X.-M., and Buries, S. 1996, Nature, 381, 207 
Viegas, S.M. & Clegg, R. 1994, MNRAS, 271, 993 



10 



Figure Captions 



Figure 1: The effect of temperature fluctuations on the hehum and oxygen abun- 
dances for our fiducial data set. The left tick mark is for no temperature 
fiuctuations (AT = 0) and the right tick mark is for AT = 2000K. 

Figure 2: The temperatures of the Hii regions in our fiducial data set as a function 
of the oxygen abundance (evaluated for AT = 0). 
Figure 3a: The distribution of offsets in the primordial mass fraction of ^He (inferred 
from a linear fit to Y versus 0/H) due to temperature fiuctuations for 
AT^„^ = 2000K (shaded histogram) and AT^„^ = 4000K. The results 
here are for the fuU data set (0/H < 1.5x10-^). 
Figure 3b: The corresponding distribution of the ratios of the slopes of the Y versus 
Z (Z 20 (0/H)) relations for ATmax = 2000K (shaded histogram) and 
AT„(a^ = 4000K compared to the fiducial slope (for AT = 0). 

Figure 4: The distribution of offsets in the primordial mass fraction of "^Hc evalu- 
ated for the full data set with AT^ax = 4000K for three different tem- 
perature probability distributions (from top to bottom: ffat, linearly in- 
creasing, linearly decreasing). 

Figure 5: In the three panels the solid curves show the variation of the mean values 
of the offsets in Yp (top panel), the ratio of slopes (middle panel) and 
the dispersion around the best fit linear Y vs. Z relation (bottom panel) 
as a function of AT^ax- The dotted curves show the 95%CL ranges. 
Figure 6a: Similar to Figure 3a, but for the low-metallicity subset of Hii regions 
(0/H < 0.9x10--^). 

Figure 6b: Similar to Figure 3b, but for the low-metaUicity subset of H ii regions. 
Figure 7: Similar to Figure 5, but for the low-metaUicity subset of H ii regions (0/H 
< 0.9x10-^). 

Figure 8: The distribution of offsets in the means of the ten lowest- Y Hii region 
helium abundances for AT^a^; = 2000K (shaded histogram) and AT^a^. 
= 4000K. 

Figure 9: Similar to Figures 5 & 7 for the means of the ten lowest- Y Hii region 
helium abundances as a function of ATm„^. 
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